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The dynamics of an initially excited two-level atom in a lossy cavity is studied by using the 
quantum trajectory method. Unwanted losses are included, such as photon absorption and scattering 
by the cavity mirrors and spontaneous emission of the atom. Based on the obtained analytical 
solutions, it is shown that the shape of the extracted spatiotemporal radiation mode sensitively 
depends on the atom-field interaction. In the case of a short-term atom-field interaction we show 
how different pulse shapes for the field extracted from the cavity can be controlled by the interaction 
time. 

PACS numbers: 42.50.Pq, 37.30.+i, 42.50.Lc 



I. INTRODUCTION 

A single atom interacting with a quantized radiation- 
field mode in a high-Q optical cavity plays an important 
role in quantum optics not only due to its conceptual rel- 
evance, but also because it appears as a basic element in 
various schemes, such as in the field of quantum informa- 
tion science (for a review see, e.g., Refs. [H, H, d, 3 ) • Cav- 
ity quantum electrodynamics (QED) has been used for 
the generation and processing of nonclassical radiation, 
as, for example, in the siiigle-atom maser 0, 0, H, [§| 
or in the optical domain 'lO'|. The quantum control of 
single-photon emission from an atom in a cavity for gen- 
erating single-photon Fock states on demand has been 
realized Jd|, and single-photon Fock state generation of 
high efficiency has been a key requirement in various 
applications such as quantum cryptography (l^ . IT3| or 
quantum networking for distribution and processing of 
quantum information [3, [l^. Recently, single-photon 
sources operating on the basis of adiabatic passage with 
just one atom trapped in a high-Q optical cavity have 
been realized [H, [Tj, [H, [l^ . In this way, the generation 
of single photons of known circular polarization has been 
possible [201 . Moreover, the adjustment of the spatiotcm- 
poral profile of single-photon pulses has been achieved 

Mill. 

In view of the widespread applications of cavity- 
assisted single-photon sources, it is of great importance 
to carefully study the quantum state of the field escap- 
ing from a cavity. Let us consider the simplest case of 
a two-level atom that near-resonantly interacts with a 
narrow-band cavity-field mode. On a time scale that is 
sufficiently short compared to the inverse bandwidth of 
the mode, the radiative and non-radiative cavity losses 
may be disregarded, and the atom-field dynamics can be 
described by the familiar Jaynes-Cummings model [23j . 
Clearly, for longer times, the atom-cavity system can no 
longer be regarded as being a closed system, and the 
losses must be taken into account. Since the wanted out- 
going field represents, from the point of view of the atom- 



cavity system, radiative losses, the study of the input- 
output problem necessarily requires inclusion in the the- 
ory of the effects of losses. Such a system, consisting of a 
two-level atom interacting with a single mode of a lossy 
cavity, has been widely considered in the past decades. 
Some of the initial theoretical works treating the effects 
of losses on the Jaynes-Cummings dynamics can be found 



in Refs. |2J, l25|, [261, |2] 
on this topic see, e.g. 




^12 
Refs. 



. For a review 
Anyway, a detailed 
characterization of the cavity output field in such a sys- 
tem still presents some open questions of significant in- 
terest, as, for example, the control of the pulse shape of 
the emitted photon. 

There are primarily two approaches to this problem, 
which are based on either quantum field theory or quan- 
tum noise theory. In quantum field theory, the system is 
commonly described on the basis of Maxwell's equations 
as used in macroscopic QED (ssl [36j . It has been shown 
that an approximate description of the fields inside and 
outside a cavity can be formulated in terms of quantum 
Langevin equations and input-output relations [37, 3§|. 
Macroscopic QED can also be used to study effects of 
unwanted losses, such as scattering and absorption losses 
caused by the cavity mirrors [H, [4l[ . More recently, 
the photon emission by an excited atom in a cavity has 
been analyzed by the method of macroscopic QED [i^] ■ 
By using a source-quantity representation of the electro- 
magnetic field, the properties of the outgoing field are in- 
vestigated. In such an approach the field inside and out- 
side the cavity is combined in a unique radiation mode, 
without regarding the fields inside and outside the cavity 
as representing independent degrees of freedom. 

Conversely, in quantum noise theory the fields inside 
and outside a cavity are regarded as representing inde- 
pendent degrees of freedom [H, [3, E^]- Accordingly, 
such a theory is based on discrete and continuous mode 
expansions of the fields inside and outside the cavity, re- 
spectively. Thus the operators of the intracavity and ex- 
ternal fields are regarded as commuting quantities. The 
continuum of the external modes is regarded as play- 
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ing the role of a dissipative system. Its effect on the 
dynamics of the intracavity modes can be treated by 
quantum Langevin equations, or, alternatively, by mas- 
ter equations [i^, [43, H^. For obtaining their solu- 
tion one can apply, for example, the quantum trajectory 
method [49l ISOL ISlj. 

In the present paper we consider a two-level atom in- 
teracting with a lossy cavity, within the framework of 
quantum noise theory, giving particular emphasis to the 
derivation of the pulse shape of the emitted photon. The 
dynamical evolution of the open quantum system under 
study is described by a master equation for the reduced 
density operator of the atom-cavity system. The effects 
of unwanted losses, such as spontaneous emissions of the 
two-level atom out the side of the cavity and photon ab- 
sorption and scattering by the cavity mirrors, are also 
taken into account. For an initially excited two-level 
atom in an empty cavity, we solve the master equation 
analytically by using the quantum trajectory method. In 
order to characterize the cavity output field, described 
by a single-photon spatiotemporal mode, we connect the 
probability to measure a photon in this mode with the 
photodetection probability as given by the quantum tra- 
jectory theory. This allows us to derive the shape of 
the mode of the extracted cavity field, which shows a 
clear mapping of the intracavity-field dynamics onto the 
mode of the output field. Moreover, the probability of 
the outgoing mode to carry a single-photon Fock state 
is calculated. After considering the case of a continu- 
ing atom-field interaction, we also analyze the short-term 
atom-field interaction. It is shown that, by changing the 
interaction time of the atom with the cavity, the mode 
structure of the outgoing field results in different pulse 
shapes of the single-photon wave packet. This opens pos- 
sibilities to control the shape of the pulse. 

The paper is organized as follows. In Sec.[n]the master 
equation describing the dynamics of the atom-cavity sys- 
tem is introduced, and the problem is solved analytically 
by using the quantum trajectory method. In Sec. Illli 
giving a description of the single-photon wave packet in 
terms of spatiotemporal mode functions, the shape of 
the mode of the extracted cavity field is obtained. In 
Sec. lIVI we analyze short-term atom- field interactions and 
obtain different shapes of the mode of the extracted field. 
A summary and some concluding remarks are given in 
Sec.lVl 



II. DAMPED ATOM-FIELD DYNAMICS 

In this section we analyze the dynamics of the sys- 
tem under scrutiny starting from a master equation and 
solving it by using the quantum trajectory theory. We 
consider a single two-level atomic transition of frequency 
uja coupled to a cavity mode of frequency uJc- The cav- 
ity mode is detuned by A from the two-level atomic 
transition frequency, = oJc + A, and is damped by 
losses through the partially transmitting cavity mirrors. 





FIG. 1: The cavity mode of frequency Uc is detuned by A 
from the two- level atomic transition frequency Li^a = i^c -I- A. 
Ki and K,2 are the photon escape rate of the cavity and the 
cavity mirrors' absorption and scattering rate, respectively. V 
is the dipole relaxation rate from level [2) to 



cf. Fig. [TJ In addition to the wanted outcoupling of the 
field, the photon can be spontaneously emitted out the 
side of the cavity into modes other than the one which 
is preferentially coupled to the resonator. Moreover, the 
photon may be absorbed or scattered by the cavity mir- 
rors. 

Treating the dissipation due to the cavity losses in a 
standard way [i^, [43, [11| , the dynamical evolution of the 
reduced density operator p[t) of the atom and the cavity 
field is described by the following master equation 



dp{t) _ 1 
dt ih 

r 

+ 2 



H, pit)] +J2 y [2a/5(i)at-atap(t)-p(t)ata] 



(1) 



2ii2/5(t)i21 - A22/5(i) - p{t)A22 



Here ki and K2 are the photon escape rate of the cavity 
and the cavity mirrors' absorption and scattering rate, re- 
spectively. We denote by F the spontaneous emission rate 
of the two-level atom. The Hamiltonian that describes 
the atom-cavity interaction is given, in the rotating-wave 
approximation, by 



H = 



hg {d 



A 



21 



12 



22 



(2) 



where d and are annihilation and creation opera- 
tors for the cavity field, respectively, and Aij = 

— 1,2), where |1) and |2) are the two atomic en- 
ergy eigenstates. Moreover, g is the atom-cavity cou- 
pling constant. Here we are considering an interac- 
tion picture with respect to Hq — hAA22, where Hq = 
hwcd^d + {ll2)nw^{A22 - ill). 

For notational convenience, in the following we will 
identify the state \a) with the state |2,0), the atom in 
the upper level, and no photon in the cavity. The state 
I a) will be considered to be the initial state of the sys- 
tem. Moreover, we will indicate with |6) the state |1, 1), 
the atom in the lower level, and one photon in the cav- 
ity. Due to photon extraction through the cavity mirror. 
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photon absorptions, or spontaneous emissions, the quan- 
tum state of the atom-cavity system is projected into the 
state |c), that indicates the state |1,0), i.e., the atom in 
the lower level and no photon in the cavity. It follows that 
the Hilbert space that describes the atom-cavity system 
under scrutiny is, in this model, simply spanned by the 
three vectors |a), |6) and |c). 

To evaluate the time evolution of the system different 
approaches can be used, see, for example, Refs. [H, [H, 
1541 . Here we have found it quite convenient to use a 
quantum trajectory approach jiol . Isol . Isij to obtain our 
analytical solutions. In this approach the dynamical evo- 
lution of the unnormalized state vector \'ip{t)), that de- 
scribes the system at time t, is governed by a nonunitary 
Schrodinger equation with a non-Hermitian Hamiltonian. 
The evolution generated by this Schrodinger equation is 
randomly interrupted, from time to time, by the action 
of collapse, or jump, operators. 

More precisely, in our specific case, considering the sys- 
tem prepared at time to = in the state |a), to determine 
the state vector of the system at a later time t, assuming 
that no jump has occurred between time to and t, we 
have to solve the nonunitary Schrodinger equation 



H'li^noit)), 



(3) 



where H' is the non-Hermitian Hamiltonian given by 

K T - 

H' = H - ih-a) a - ih-A22 , (4) 

with H given by Eq. ([2]), and where we have defined 

(5) 



K = Kl 



K2 



If no jump has occurred between time to and t, the system 
evolves via Eq. ([3]) in the unnormalized state 



iV^noW) -a(t)|a)+/3(t)|6) 



(6) 



In this case the conditioned density operator for the 
atom-cavity system is given by 



/5no(t) 



|^no(t))(^no(t)| 
(^no(t)|^„o(t)) 



(7) 



Here we have used the word conditioned to stress the 
fact that this is the density operator at time t one ob- 
tains conditioned to the fact that no jump has occurred 
between time to and t. 

The evolution governed by the nonunitary Schrodinger 
equation ([3]) is randomly interrupted by three kinds of 
jumps, Ji, J2 and Js given by 



Ji = \/^a [i = 1, 2) . 

is = Vf il2. 



(8) 
(9) 



The jump operators Ji and J2 are related to a photon 
extracted from the cavity and a photon absorbed or scat- 
tered by the mirrors, respectively. The jump operator 



Js is related to a photon spontaneously emitted by the 
atom. If a jump has occurred at time tj, tj G (to,t], the 
wavevector is found collapsed in the state |c) due to the 
action of one of the jump operators 

J^\i>no{h))=^J^^a\i^no{ti)) ^\c) (i = l,2), (10) 
Js |l/i„o(t,j))=Vril2|V5no(tj)) ^ |c) . (11) 

It is clear that in the problem under study we can 
have only one jump. Once the system collapses in the 
state |c) the nonunitary Schrodinger equation ([3]) simply 
keeps it there forever. In this case the conditioned density 
operator at time t is given by 



Pycs(t) = |c)(c| 



(12) 



where we indicate with "yes" the fact that a jump has 
occurred. 

According to the quantum trajectory method, the den- 
sity operator p(t) is obtained by performing an ensemble 
average over the different conditioned density operators 
at time t. In the present case, starting at time to with 
the density operator po — \a){a\, the ensemble average is 
performed over the two possible realizations (histories) 
"yes" and "no": 



/5(t) = Pno(t)/5no(t) + Pycs(t)/5ycs(t) 



(13) 



Here Pno(t) and pycs(t) are the probability that between 
the initial time to and time t no jump and one jump has 
occurred, respectively. Of course, Pno(t) + J5ycs(t) — 1. 
The density operator given by Eq. (fTS]) tells us that the 
system at time t is in a statistical mixture: either no 
photon has escaped from the cavity or one (and only 
one) photon has escaped. 

To evaluate Pno(t) we use the method of the delay func- 
tion [50]. This method tells us that the probability pno(t) 
is given by the square of the norm of the unnormalized 
state vector: 

P„o(t) =11 \4>no{t))\?= (^„o(t)|V^„o(t)) = |a(t)|2-f |/3(t)|2. 

(14) 

From Eqs. (fT3|) and (fT4| one obtains for the density op- 
erator /5(t) the expression 

p{t) = \a{t)\^\a){a\ + \(i{t)?\h){h\+a{t)P*{t)\a){b\ 
+ a*{t)m\b){a\ + \l{t)?\c){c\, (15) 



where we have defined 



|7(t)P^Pyes(t) = l-[Kt)P + |/3(t)| = 



(16) 



The physical meaning of |a(t)p, |/3(t)p and |7(t)P is 
clear. They represent the probability that at time t the 
system can be found either in |a), \b), or |c). Moreover, 
from the master equation ([1]) , together with Eq. pS]) , one 
obtains 



d\l{t)? 
dt 



= Tr 



dm 

dt 



c){c\ 
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To better understand the meaning of Eq. (|T7|) , it is use- 
ful to do the following consideration. According to the 
quantum trajectory theory, the probability for a jump, cf. 
Eqs. ^ and to occur in the time interval {t,t + dt] 
is given by (i = 1, 2) 

Pi{t) = {jjj,)tdt = KiTr [p{t)d^a] dt ^ K,\(3{t)\^dt (18) 
and 

p,{t) = {,P,J,)tdt = TTr [p(t)i22]dt = r|a(i)pdt. (19) 

Of course, the increment in the time interval dt for pycs{t) 
is equal to pi{t) +p2{t) +ps{t) , so that we can write, using 
Eqs. (0, (IinD, and 

d\-fit)\^ ^ dpy,,{t) ^ K\(3it)\^dt + T\ait)\^dt , (20) 

that is again Eq. (flT)) . The physical meaning of this re- 
lation is quite clear. When the system is in \b), i.e. with 
probability |/3(t)p, we can have an emission of a pho- 
ton from the cavity or an absorption or scattering by the 
cavity mirrors (controlled by the parameter k). When 
the system is in \a), i.e. with probability |a(i)p, we can 
have a photon spontaneously emitted by the atom (con- 
trolled by the parameter F). The related jumps operators 
project the system into |c), hence producing an incre- 
ment of |7(i)p. Moreover, by integrating equation (I^Dl) 
one gets 

Pycs 

(t) = \'y{t)\ = pcxt(i) +Pabs(i) + Pspo (t), (21) 
where we have defined 



and 



Petit) ^Kl / dt'\P{t')\\ 
JQ 

Pabs(t) =K2 / dt'\P{t')W 

Jo 



Pspo(t) = r / dt'\a{t 



(22) 
(23) 

(24) 



The function Poxt(i) represents the probability that a 
photon is extracted from the cavity in the time inter- 
val [0,t], and Pa,hs{t) the probability that a photon is 
absorbed or scattered by the mirrors in the same time 
interval. Finally, Pspo{t) represents the probability that a 
spontaneous emission has occurred in time interval [0, t]. 

Note that from Eqs. (HH), and ^ it follows that 
Pcxt{t), Pahs{t), and Pspo(i) have to be monotonically in- 
creasing functions: the longer one waits, the larger is the 
probability that a photon has leaked out of the cavity, or 
is absorbed or scattered by the mirrors, or a spontaneous 
emission has occurred. Moreover, if we wait long enough 
a photon is certain to be emitted in one of the three ways, 
so that limt^oo l7(i)P = 1- In this case Pcxt(i) does not 
reach asymptotically the value 1, due to the presence of 
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FIG. 2: The probabilities \a{t)\^ (dashed line), \f3{t)\^ (dotted 
line), pext(t) (full line), and Pabs(i), Pspo{t) (dot-dashed lines) 
are shown for 2g/K — 10, A/k — 0.1, ki/k — 0.9, and F/k = 
0.5. 



spontaneous emissions and mirror absorption or scatter- 
ing. If we take the limit f ^ oo of Eq. pT|) , we get 



Pcxt(oo) + Pabs(oo) -I- Pspo(oo) = 1 , 



(25) 



In order to determine a{t) and (3{t) we have to solve the 
nonunitary Schrodinger equation, cf. Eqs. Q and ([4]). 
This brings us to consider the following linear system of 
differential equations, 

.■r\ 



r a{t) = -i (A - i^) a{t) - igl3{t) 
{ (3{t) = ~iga{t) - f /3(t) . 



(26) 



For the initial conditions a(0) = 1 and /3(0) = 0, and 
defining 



IK 



A- 



(27) 



we can write the solutions as 

a{t) 



K/2-i(A-iT/2) , fflt 

smh — 

Q V 2 



cosh ( — 
2 



-[(K+r)/4-(-jA/2]t 



P{t) = -^ sinh f ^ ) e-[(-+n/4+^A/2]t (28) 



nt 



Using the solutions given in Eq. (pS)) . one can plot the 
probabilities to find at time t the system in \a) or in \b), 
i.e. |Q!(t)p and |/3(t)p, respectively, as well as Pcxt{t), 
Pabs(t) and ftpo(i), as given by Eqs. (I121)-(E11)- In Fig. [2] 
we show these functions for the parameters 2g/K — 10, 
A/k = 0.1, ki/k = 0.9, and for the realistic choice of 
F/k = 0.5, cf. Ref. H. Note that for t oo, \a{t)\^ 

0, |/3(0P ^ 0, and pext(i) +Pabs(0 +Pspo{t) = |7(<)l' ^ 

1, as it is expected from Eqs. ([^5]) and fTB]) . 

Let us consider A = and g ^ k, F. From Eq. 
one immediately obtains 

-(K+r)t/2 



\m\' 



cos^ {gt)e 
sin2(5i)e-(''+r)*/2 



(29) 
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In this case the system undergoes damped Rabi oscil- 
lations between \a) and \b) with frequency 2g. From 
Eq. ((291) one easily obtains, using Eqs. (|22l) - ([24ll . 



pUt)^^il dmtT = ^ (i_e-('^+r)*/2) , (30) 
P.Ut) = n,J^dmt'r = -!^ (l-e-('^+n*/2) ,(31) 



and 



1-e 



-(K+r)t/2 



(32) 



Pspo(t) = r/ dt'\a{t')\^ = ^ 

Jo K -\- L 

From Eq. (|2ip we then get 

Pyes(^) = |7(^)l' = l-e-(^+^)'/^ (33) 
that shows a simple exponential behavior. For t—KX) we 

have Pcxt(t)^Kl/(K + r), Pabs(0— ^«2/(K + r), Pspo(t)^ 

r/(K + r), and |7(<)|2^1. 



and from Eq. ((361) . it follows that 

'b{t),b\t')\ ^S{t-t'). 



(38) 



Given that the photon is in the mode ^i, i.e. it is de- 
scribed by the normalized function £,i(t), according to 



(39) 



it is possible to construct a complete orthonormal set 
{^j{t)} of functions where 



and 



(40) 



(41) 



III. SINGLE-PHOTON WAVE PACKET 

The analysis performed here, using a quantum trajec- 
tory approach, is implicitly based on an unraveling of 
the master equation ([1]) for the case of direct photoelec- 
tric detection of the field emitted from the cavity [5l| . 
In experiments one uses a large number of photodetec- 
tions to recover the properties of the electromagnetic 
field, also in the case of single-photon sources, see, e.g., 
Refs. [Ill, [22|. Properties like wave packet duration or 
wave packet bandwidth are analyzed with an ensemble 
of photons and cannot be determined from a measure- 
ment on just a single photon. In this respect it is impor- 
tant to carefully describe the arrival of photons at the 
photodetector. 

For this purpose it is convenient, following the ap- 
proach of [iaHj,!!!!, to choose spatiotemporal modes for 
characterizing the single-photon wave packet. A stream 
of single photons emitted one after the other can be de- 
scribed by a state vector 1 1^ . ) , where 



4lo>. 



(34) 



Here c|. is the creation operator for photons of spatiotem- 
poral mode £,i{t) defined as 



where ^i{t) — for t < 0, and 6^(t) is given by 
1 



dujb\uj)e'' 



(35) 



(36) 



that is, the Fourier transform of the operator ¥{ijj), the 
creation operator of quanta of a monochromatic wave of 
frequency a; in free space. From the relation 



b{uj),b\uj') \ =5{uj-J). 



(37) 



From Eqs. ([38]) and (|40l) it is immediate to show that 



(42) 



so that the operators defined by Eq. (|35p using the com- 
plete set of orthonormal functions represent a set 
of independent bosons, and c^. can be used to construct 
number states in the usual way. 



1 



|0). 



(43) 



The spatiotemporal mode function ^^(i) is composed 
of an amplitude envelope €i{t) and a phase 4'i['^)i 

^,(i) = e,(t)e'*'(*). (44) 
According to Eq. ((55)) . the normalization reads as 



dtm = i. 



(45) 



If we now define the flux operator in units of photons per 
unit time, 



f{t) = b\t)b{t), 



(46) 



using the inverse relation of Eq. ((55|) . b{t) — 
we can write Eq. ((46)) as 



(47) 



When no extra losses, such as spontaneous emissions out 
the side of the cavity or mirrors' absorption or scattering, 
are considered, the density operator of the cavity output 
field for a photon in the mode is, in the Heisenberg 
picture, pout = llfJil^iL cf. Ref. (58||. When extra losses 



6 



are included, the density operator of the cavity output 
field is given by the statistical mixture 

Pout =Pext(00)|l5J(l,..| + [1 -Pext(oo)] |0)(0| . (48) 

Note that 1 — Pcxt(oo) = Pabs(oo) +Pspo(oo)i consistently 
with the fact that the zero-field contribution is related 
to the spontaneous emissions out the side of the cavity 
or to mirrors' absorption or scattering. The probability 
density distribution of measuring the photon at a given 
time t is then 



P^.(t) = Tr pout/W =Pext(oo)e^(i) 



(49) 



Of course, integrating Eq. (|^ . and using Eq. we 
have 

-Ptot = / dtP^.{t) ^ pc^t{oo) / ej(t) = pcxt(oo) , 
Jo Jo 

(50) 

as it is expected. 

Let us consider a photon in the mode fi, whose ampli- 
tude envelope ei{t) does not change significantly in the 
detection time resolution T. The response probability of 
the detector within a time interval [t — T/2, t + T/2] is 
then, using Eq. given by 



rt+T/l 
Jt-T/2 



ioxt(oo)e,^(t) T. 



(51) 



In the case of a detector of quantum efficiency 77, Eq. ([5T|) 
becomes 



-Pd(0 ='7Pcxt(oo)e2(t)T. 



(52) 



In a usual experiment a large number of photodetection 
events are accumulated to obtain P£,(t), and from these 
measurements one gets ei{t). This consideration is im- 
portant because it tells us how one can obtain the am- 
plitude envelope of the mode function within a quantum 
trajectory formalism. The probability to measure be- 
tween time t — T/2 and t + T/2 a "click" at the detector 
is equal to the probability to have a jump Ji in the same 
time interval, so that using Eq. (fT8|) . we get, in the case 
of a detector of efficiency 77, 

PD(t) - r/«iTr [/5(t)ata] T = ryKi|/3(i)|2r . (53) 

Comparing this with Eq. (|52p we obtain 



Poxt(oo) 



1/3(01 



(54) 



Let us check if Eq. (gSl) is still fulfilled. Using Eq. 
we obtain 



/ dteUt)^ . . 

Jo Pcxt(oo) t 



lim Kl / dt'\l3{t')\^ = 1. (55) 




To obtain Eq. ([54|) we have assumed that the photode- 
tector was positioned just outside the cavity, at 2; = 0, 




FIG. 3: Plot of the amplitude envelope ei{z, t) \/pcxt (00) /ki 
for the spatiotemporal mode of the cavity output field, with 
parameters 2g/K = 10, A/n — 0.1, hi/k = 0.9, F/k = 0.5, 
and Kt = 7. 



SO that this equation gives the amplitude envelope of the 
spatiotemporal mode function at z = 0. Of course we 
could imagine to position the detector in an arbitrary 
position z (with z > 0). In this case, because outside 
the cavity the field is propagating at the speed of light 
c, the amplitude envelope ti{z,t) is related to the one 
at position z = via the retarded time t ~ z/c. More 
precisely. 



z(0,t- 



Pcxt(oo) 



|/3(i-f)| i>f>0 







t < 



(56) 

Note that |/3(i)P represents the probability to find at 
time t a photon inside the cavity. In this respect, Eq. ([56]) 
shows that the intracavity field dynamics determines the 
structure of the spatiotemporal mode of the output field. 
Moreover, in order to clarify the connection between 
ei{z, t) and the probability to emit a photon, it is useful to 
consider the following equation, obtained from Eqs. (I5(: 
and dig): 



z/c 



dt'eUz,t') = P-'^'-'('^ 

Poxt(oo) 



{t > z/c) , (57) 



where pcxt (t) represents the probability that a photon has 
leaked out of the cavity in the time interval [0,t]. Using 
Eq. dnSl) and the solution given by Eq. ^E^, in Fig. [3] 
we plot the amplitude envelope ei(z, t)-\/j3ext(oo)/Ki for 
the spatiotemporal mode of the cavity output field for 
the case where the parameters g and k are chosen as 
2g/K = 10, A/k = 0.1, T/k = 0.5, and Kt = 7. It 
is clearly seen that the intracavity dynamics strongly 
modulates the mode structure of the photon wave packet 
propagating outside the cavity. 
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IV. TIME CONTROL OF THE PHOTON WAVE 
PACKET 

Let us now consider the case where the interaction of 
the atom with the cavity-assisted field has a hmited du- 
ration, so that it effectively terminates at time r. To 
analyze this situation we have obviously to split the dy- 
namical evolution of the system in two distinct time in- 
tervals, one interval from the initial time io to the 
time T, and the second interval for times t greater than 
T. For times t such that < t < r, the evolution is still 
described by the one previously analyzed in Sec. |TT1 with 
solutions given by Eq. ([^5]) . For times t with t > t, when 
the interaction of the atom with the cavity is set to zero, 
the Hilbert space of the system (cavity field) reduces to 
that spanned by the two Fock-state vectors |1) and |0) of 
the cavity field. At time t the cavity field, obtained by 
tracing over the atomic states in Eq. ([T5)) . is described 
by the following density operator: 

p(r)=[l-|/3(r)p]|0)(0|-|-|/3(r)ni)(l|, (58) 



where /3(r) is given by Eq. 

To analyze the dynamical evolution, for t > r, of this 
initial state, we follow the procedure given by the quan- 
tum trajectory theory when the initial state is not a pure 
state, but a statistical mixture [iiliillH. For the part 
related to [l - |/3(r)p] |0)(0|, one has to start with the 
state vector \iP{t)) = |0). Because in this case the cavity 
is already empty, the evolution simply leaves the cavity 
in its vacuum state also at later times, so that we simply 
have 



Pi(t) = |0)(0| 



(59) 



For the part related to |/3(r)p|l)(l|, one has instead 
to start with the state vector \iIj{t)) = Before 
a collapse occurs, the evolution of the unnormalized 
state |-0no(i)) = P{t)\l) is described by the nonunitary 
Schrodinger equation 

«^^|^„o(t)) =-ih'^a^d\iJr,oit)). (60) 

Its solution is /3(t) = e-''(*-^)/2_ if a jump J,, cf. Eq. dH]), 
has occurred at time tj, tj € (r, t] , the wave vector col- 
lapses in the state |0), 



i»|^no(ij)) = ^^a[P{tJ)\l)] 



(61) 



This tells us that at time t > t the density operator for 
the part related to the initial state |^(t)) = |1) is, using 
the method of the delay function [53|, given by 



P2{t)^Pnom){l\ + [l-Pno{t)] |0)(0| , (62) 



where 



Pno{t)=\\ \i^n,{t))\\'^\m\' 



-K(t-T) 
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FIG. 4: The probabilities Pin(i) (dotted line) and Poxt(i) (full 
line) are shown for r = ■k/\Q.\, 2g/hi — 10, A/k — 0.1, = 
0.9, and T/k = 0.5. 



|0)(0|. (64) 



Substituting this in Eq. (|62|) one gets 
P2(t) = e-«(*-^)|l)(l|+ [l-e-''(*-^) 



The evolution for the initial density operator given in 
Eq. (|58p is, according to the quantum trajectory theory, 
given by 

m = [i-mr)\']pi{t) + mr)\'p2it), m 

so that we obtain, using Eqs. ([55)1 and ([M]) . the result 

m = |/3(r)pe-'^(*-)|l)(l| 



(66) 



This is the density operator for the intracavity field for 
t > T. Combining this result with the one for the time 
interval [0,t], cf. Eq. (fT5|) . we can write the probability 
to find a photon inside the cavity at an arbitrary time t 
as 

Pi„(i) = e(r - t)\m\' + 0(t - r)|/3(T)|2e-«(*--), (67) 

and e(t) 



where |/3(i)p and |/3(r)p are given by Eq. 
is the unit-step function. 

Let us now consider the probability pcxt (t) that a pho- 
ton is extracted from the cavity in the time interval [0,t\. 
For t < T this probability is equal to pcxt(^), cf- Eq. (|22p . 
with (3{t) given by Eq. P5|) . For t > t we have a sum 
of two contributions, a first one up to time r, given by 
PcxtiT), and a second one for the time interval [r, t]. This 
second contribution is given, using Eqs. (HH]) and by 



dt'{jl,h)t' = Ki I dt'\ 



(^)|2g-K(t'-r) 



Kl 



1 - e 



-K(t-r) 



(68) 



We can now combine these results and write, for an ar- 
bitrary time t, the probability Pcxt{t) as 



(63) Pext(t) = e(T-t)pext(t) + e(t-r)[pcxt(T)+p:,t(t)]. (69) 
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From this equation one gets that, ior t oo, the ex- 
traction probabihty is equal to Pext(oo) = Pcxt(T) + 
{ki/k)\(3{t)\^ . This relation can be rewritten, using 
Eqs. dill) and dUD, as 



Poxt(oo) = 1 - |Q;(r)p -_Pspo(t) -_Pabs(oo) 



(70) 



where Pabs(oo) = Pabs(T) + (k2/k)|/3(t)|2, and Pabs(T), 
Pspo(t) are given by Eq. (^5)1 and Eq. respectively. 
Eq. l|70|) shows that the extraction probability is not, 
in general, asymptotically reaching the value one. This 
reflects the fact that at time r the atom can be, with 
probability |a(T)p, in its excited state. If the interaction 
is set to zero when the atom is in its excited state, then, 
obviously, no photon extraction can anymore occur from 
the empty cavity. Moreover, also the contribution due to 
spontaneous emissions up to time r, pspo('''), and the to- 
tal absorption probability, Pabs(oo), remove photons from 
the extracted output channel. Note that from Eq. ([70]) 
one obtains, for t oo, Eq. (pS)) . as it is expected, being 
l«(r)P-0. 

Let us now analyze the dynamical evolution of the sys- 
tem for the following two cases. In the first case we con- 
sider that the atom leaves the cavity at time r equal to 
the first half Rabi cycle, i.e. at r = 7r/|17|, cf. Eq. (p8| . 
In the second case we consider that the atom leaves the 
cavity at time r equal to the first Rabi cycle, i.e. at 
T = 2n/\n\. For t = 7r/|rj|, using Eqs. §7}i and dM]), we 
plot in Fig. |4] the behavior of the probabilities Pin{t) and 
^ext(i); respectively. Note that for the parameters used 
in this case we have |a(r)p ~ 0. For r = 27r/|ri|, the 
behavior of the probabilities Pin{t) and Pext(i) is plotted 
in Fig. [5l using Eqs. (|67| and ([69)l . Because in this case 
|/3(t)P ~ 0, one sees from Eq. (p7)) that the probability 
to find a photon inside the cavity at time t > r is con- 
stant and, approximately, equal to zero. From Eq. ([M]) 
the probability Pcxt(i) has, for t > t, the constant value 
Pcxt(T). Moreover, in this case |a(T)p is not negligible 
because the interaction was switched off when the atom 
had a significant probability to be found in its excited 



state. 

We consider now the problem to determine the single- 
photon pulse shape in the case of a short-term atom- 
field interaction. Because for t > t the probability to 
measure between time t — T/2 and t + T/2 a "click" at 
a photodetector of efficiency 77, positioned at 2: = 0, is 
given by Pvit) = 7]Ki\P{t)\'^ e~'^^-*~^^T , we can write the 
amplitude envelope of the spatiotemporal mode function 



/ pext(oo) ^ 



,{t) = e{T-t)\/3{t)\+Qit~r)\l3{T)\e 



where p{t) and /3(t) are given by Eq. 



(71) 

, and Poxt(oo) 

by Eq. ([70|) . The factor \/pcxt{oo) / ki is needed in or- 
der for the spatiotemporal mode function to be properly 
normalized accordingly to Eq. (|45p . If we now define the 
retarded time t,- = t — z/c, we can generalize Eq. (j7ip as 



/ pext(oo) 



+ e(i,-T)|/3(T)|e 



-f (*r-r) 



, (72) 



for t > z/c > 0, and, obviously, ei{z,t) = for t < z/c. 

Using Eq. ([7^ one can obtain, for example, the ampli- 
tude envelope of the spatiotemporal mode for the two 
cases above considered, i.e. for r = t^/I^I and for 
T = 27r/|f2|, respectively. In Fig. [H] we consider the case 
where the interaction time r is chosen as r = 7r/|f2|, and 
Kt = 7. After |/3(i)p, the probability to find a photon 
inside the cavity reaches its maximum values at time r, 
the interaction with the atom is set to zero, and one 
observes an exponential decay regulated by the photon 
escape and absorption rate k. This behavior is mapped 
in the amplitude envelope shape of the extracted spa- 
tiotemporal mode, as can be clearly seen in Fig. [Sj If 
we consider now the case where the interaction time r 
is chosen as r = 27r/|f2|, using Eq. ([7^ we obtain, for 
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FIG. 7: Plot of the amplitude envelope ei (2, t) \/pcxt (00) / ki 
for the spatiotemporal mode of the cavity output field, for 
r = 27r/|r2|, = 7, = 10 and A/k = 0.1, k.i/k = 0.9, 

and T/k = 0.5. 




17 18 19 20 



FIG. 8: Plot of the amplitude envelope ei{z, t) •\/pext(oo) /ki 
for the spatiotemporal mode of the cavity output field in the 
region 17 < kz/c < 20, for kt = 2.2, Kt = 20, 2g/K = 10, 
A/k = 0.1, ki/k = 0.9, and for T = (fuU line), or F/k = 0.5 
(dashed line). 



Kt = 7, the function plotted in Fig. [71 Here, |/3(t)p, after 
reaching its maximum value, returns practically to zero, 
for T = 27r/|r2|. At this point the interaction with the 
atom is switched off, so that the cavity is, practically, 
left with no photon inside, so that no photon can be ex- 
tracted at later times. This dynamics is clearly mapped 
in the amplitude envelope shape of the spatiotemporal 
mode, producing a short pulse, of length 27rc/|ri|, as can 
be seen in Fig. [T] 

Finally, we consider the amplitude envelope for the 
spatiotemporal mode of the cavity output field for 
an interaction time r arbitrary chosen so that kt = 
2.2. In Fig. [HI using Eq. ((721), we plot the function 
ei{z, t) pcxt{oo) / Ki in the region 17 < kz/c < 20. In 
this plot the amplitude envelope shows a behavior that 
is intermediate between the one depicted in Fig. [5| and 
the one in Fig. [3] After approximately three and a half 
Rabi cycles, the interaction is switched off, and the cav- 
ity field simply decays with an exponential behavior reg- 
ulated by K. To clearly see the effects of spontaneous 
emissions we also plot in Fig. [8] the behavior of the func- 
tion ei{z, t) y^pcxti'X)) / Ki for the case F = 0. We see that 
the presence of spontaneous emissions is not negligible. 

In order to realize the time control analyzed here, let 
us consider a neutral atom [3, [l^ or a trapped ion 
in an optical cavity. With the use of an external laser 
pulse it is possible to excite the atom to an auxiliary elec- 
tronic state, to decouple the atom from the cavity mode. 
In this way the interaction time between the atom and 
the cavity can be regulated. Moreover, for a continuous 
time-dependent control of the interaction, a pulsed Ra- 
man coupling could be useful, or the atom could be tuned 
out of resonance by external electric or magnetic fields. 
A single-photon wave packet with a defined pulse shape 
may be used for a bidirectional atom-photon interface 
in a quantum network [T^ . This relies on the coher- 
ent interaction between the atom and the cavity field, 
provided that the effective coupling exceeds the atomic 



decay rates [soj . 



V. SUMMARY AND CONCLUSIONS 

The dynamics of an atom-cavity system, which consists 
of an initially excited two-level atom in a lossy cavity, has 
been analyzed. The open quantum system under study 
has been described by means of a master equation. By 
using the quantum trajectory method, we have derived 
analytical solutions for the dynamics of the system. The 
effects of unwanted losses, such as absorption and scat- 
tering by the cavity mirrors and spontaneous emission 
of the atom into field modes out the side of the cavity, 
have also been considered. After giving a description of 
the single-photon wave packet in terms of spatiotempo- 
ral mode functions, we have connected the probability 
to measure a photon in a definite mode structure of the 
output field with the photodetection probability. In this 
way the shape of the mode of the extracted cavity field 
has been obtained. The mode shape sensitively depends 
on the atom-field interaction, showing a clear mapping 
of the intracavity field dynamics onto the output field. 
The probability of the mode to carry a one-photon Fock 
state has been calculated. We have also shown that dif- 
ferent pulse shapes of the extracted field can be generated 
by controlling the duration of the atom-field interaction 
time. 

Finally we would like to comment on a fundamental 
difference of the quantum noise approach under study 
in comparison with the treatment of the same problem 
by macroscopic QED. In the latter method one does 
not distinguish between the intracavity and the exter- 
nal fields. There exists only a unique field mode, which 
covers both the areas inside and outside the cavity. In 
the quantum noise theory on the other hand, the input- 
output coupling is introduced via the interaction of two 
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types of modes, describing the intracavity and the ex- 
ternal fields. These modes belong to different Hilbert 
spaces and are therefore commuting. In the description 
of a unified mode, as in macroscopic QED, there is no 
hint of the existence of commuting field observables that 
might describe intracavity and external fields. To over- 
come this basic difference between the two treatments, it 
seems reasonable to assume that entanglement between 
the commuting field modes in the quantum noise theory 
may replace the noncommutativity in the macroscopic 



QED in some respect. A careful study of this problem 
requires further investigations. 
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